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Abstract. 

The  present  study  examines  the  thermo-oxidative  behavior  of  high-temperature  polymer 
matrix  composite  (HTPMC)  materials.  Thermo-oxidative  aging  in  neat  resin  is 
simulated  with  a  diffusion-reaction  model  in  which  temperature,  oxygen  concentration 
and  weight  loss  effects  are  considered.  The  thermo-oxidative  behavior  of  the  composite, 
on  the  other  hand,  is  significantly  different  from  that  of  the  constituents  as  the 
composite  microstructure,  including  the  fiber/matrix  interphase/interface,  introduces 
anisotropy  in  the  diffusion  and  oxidation  behavior.  Unit  cell  analyses  are  carried  out 
using  three-dimensional  finite  element  analysis  of  repeated  volume  elements 
representing  the  fiber,  matrix  and  interphase  regions,  and  the  resin  oxidation  model. 
Parametric  studies  illustrating  the  anisotropy  in  the  oxidative  region  growth  and  the 
effect  of  fiber  and  interphase  diffusivity  on  the  oxidation  layer  growth  are  discussed.  It 
is  suggested  that  fiber-matrix  debonding  could  provide  additional  diffusion  paths  to 
explain  the  extent  of  observed  anisotropic  oxidation  growth. 
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1.  Introduction 


Polymer  matrix  composites  (PMCs)  used  in  aerospace  high  temperature 
applications,  such  as  turbine  engines  and  engine  exhaust  washed  structures,  have  limited 
life  due  to  environmental  degradation.  Predicting  the  service  life  of  composite  structures 
subjected  to  mechanical,  high  temperature,  moisture,  and  corrosive  conditions  is 
challenging  due  to  the  complex  physical,  chemical,  and  thermo-mechanical  mechanisms 
evolving  with  time  in  the  constituent  phases.  While  the  time-dependent  physical, 
chemical,  and  damage-induced  degradation  mechanisms  have  been  studied  for  some 
resin  systems,  polymer  composite  thermal  oxidation  studies  from  a  mechanistic 
perspective  are  nascent.  Notable  exceptions  to  this  are  the  recent  works  of  Colin  and 
Verdu  [1],  Colin,  et  al.  [2],  Wang  and  Chen  [3],  Pochiraju,  et  al.  [4],  Tandon,  et  al.  [5], 
and  past  work  by  Skontorp,  et  al.  [6],  Wise,  et  al.  [7],  and  McManus,  et  al.  [8],  Equally 
important  are  the  experimental  characterization  efforts  by  such  groups  as  Bowles,  et  al., 
[9];  Tsuji,  et  al.,  [10];  Abdeljaoued  [11],  and  Schoeppner,  et  al.  [12],  Current  emphasis 
is  on  the  implementation  and  extension  of  multiscale  models  to  represent  the  polymer 
behavior/properties  as  a  function  of  the  degradation  state  to  include  chemical 
degradation  kinetics,  micromechanical  models  with  time-dependent  polymer 
constitutive  relationships,  and  on  ply-level  and  laminate  level  models  to  predict 
structural  behavior.  The  behavior  of  the  composite  will  thus  be  dependent  on  its  current 
chemical,  physical,  and  mechanical  state  as  well  as  its  service  history.  This 
multidisciplinary  modeling  approach  will  provide  generalizable  analytical  and  design 
tools  for  realistic  prediction  of  performance,  durability,  and  use  life  of  HTPMCs. 
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2.  Experimental  Characterization  and  Observations 


The  attributes  of  the  oxidation  process  of  high  temperature  polymer  and  polymer 
composites  vary  dependent  on  the  chemistry  of  the  polymers  being  tested.  The  primary 
focus  of  the  discussion  in  this  work  is  on  predicting  isothermal  oxidative  aging  of 
unitape  laminates  with  particular  focus  on  the  composites  utilizing  PMR-15  high 
temperature  polymer.  Additionally,  the  reinforcement  is  limited  to  PAN-based  carbon 
fibers.  PMR-15  is  a  widely -used  addition  polyimide  with  a  maximum  service 
temperature  of  approximately  288°C.  Among  the  class  of  high  temperature  polymers  are 
the  bismaleimides,  Avimid-N,  thermosetting  polyimides  (AFR700B,  LARC  RP46)  and 
phenylethynyl-terminated  polyimides  (PETI-5,  AFR-PE-N)  resin  systems.  Each  of  these 
material  systems  has  unique  degradation  reactions,  mechanisms,  and  kinetics  particular 
to  their  chemical  structure.  Although  experimental  observations  and  predictions  of  the 
behavior  of  PMR-15  neat  resin  and  composites  are  not  necessarily  representative  of  the 
behavior  of  these  other  high  temperature  polymer  material  systems,  the  modeling 
methodology  and  procedures  for  determining  material  parameters  may  be  directly 
applicable  to  predicting  their  behavior. 

2. 1  Resin  Oxidation 

Recent  studies  [1,3,9,13,14]  document  the  growth  of  the  thermo-oxidative  layer  in 
PMR-15  resin  and  the  changes  in  elastic  moduli  and  chemical  composition  resulting 
from  isothermal  aging.  It  has  been  observed  for  PMR-15  that  the  polymer  degradation 
occurs  mainly  within  a  thin  surface  layer  that  develops  and  grows  during  thermal  aging. 
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In  our  work,  for  both  the  neat  polymer  and  composite  specimens  that  were  used  to 
monitor  the  propagation  of  oxidation,  samples  were  dry-sectioned  from  the  aged  larger 
specimen,  as  shown  in  Fig.  1,  potted  in  828-D230  epoxy  resin,  and  cured  at  room 
temperature  for  three  days.  This  sectioning  procedure  allowed  monitoring  of  four  of  the 
exposed  free-surfaces  of  the  large  specimen.  The  large  specimen  was  subsequently 
placed  back  in  the  oven  to  continue  the  aging  process.  The  diamond  blade  used  to 
section  these  samples  was  washed  with  acetone  and  wiped  clean  with  paper  towels  prior 
to  cutting  to  minimize  the  amount  of  contamination  from  the  cutting  wheel.  The 
specimens  were  then  wet-sanded  with  600-grit  sandpaper  and  distilled  water  and 
polished  using  a  0.3-micron  alumina  polishing  media.  Since  the  oxidized  layer  forms  on 
all  exposed  free-surfaces  and  propagates  to  the  interior  of  the  sample,  the  oxidized 
region  is  easily  seen  in  the  cross  section. 

Fig.  2  shows  a  photomicrograph  of  a  PMR-15  neat  resin  specimen  isothermally  aged 
in  air  under  elevated  pressure  (0.414  MPa)  at  288°C  for  a  period  of  175  hrs.  The  figure 
clearly  shows  the  oxidized  region  (much  like  a  picture  frame)  on  the  two  adjacent 
exposed  free  surfaces  of  the  specimen.  The  dark  shaped  region  between  the  outer 
oxidized  layer  and  interior  unoxidized  region  is  a  transition  region,  which  is  the  active 
“reaction”  or  “process”  zone.  These  observations  allow  easy  measuring  and 
characterization  of  the  oxidized  layer  and  the  transition  regions  as  a  function  of  aging 
temperature  and  time.  Fig.  3  shows  the  evolution  of  the  oxidized  layer  and  the  active 
“reaction”  zone  thicknesses  as  a  function  of  aging  time  for  PMR-15  resin  aged  at  288°C. 
The  thickness  of  the  oxidized  layer  is  seen  to  approach  a  plateau  value  as  the  oxidation 


4 


growth  rate  reduces  considerably  for  longer  aging  time  periods,  whereas  the  thickness 
of  the  active  “reaction”  zone  remains  nearly  constant  for  the  aging  times  considered. 

The  rate  of  oxidation  is  sensitive  to  the  partial  pressure  of  oxygen  at  the  specimen 
surface  and  acceleration  can  be  achieved  by  increasing  the  air  pressure  and  thus  the 
partial  pressure  of  oxygen  within  the  aging  chamber  (e.g.,  Ciutacu,  et  al.  [15],  Tsotsis,  et 
al.  [16]).  Recent  work  by  Tandon,  et  al.  [17]  examined  the  use  of  elevated  pressure  in 
conjunction  with  a  realistic  use  temperature  to  accelerate  the  rate  of  thermo-oxidative 
degradation  in  PMR- 1 5  resin.  It  was  observed  that  the  elevated  pressure  aging  of  the 
PMR- 1 5  neat  resin  presumably  had  a  significant  effect  on  the  rate  of  diffusion  of 
oxygen  into  the  specimen,  accelerating  the  oxidation  process  and  allowing  the  oxygen 
to  diffuse  deeper  into  the  interior  of  the  specimens.  This  results  in  greater  oxidation 
layer  thicknesses  than  are  achievable  in  ambient  air  pressure  environments  as  shown  in 
Fig.  4. 

2.2  Fiber  oxidation 

The  effects  of  service  temperatures  and  oxidation  on  the  mechanical  properties  of 
high-temperature  polymer  composites  are  primarily  manifested  in  the  polymer- 
dominated  properties,  namely,  the  transverse  properties  (perpendicular  to  the  fiber 
direction)  and  the  shear  properties.  However,  the  fiber-dominated  properties  may  be 
affected  by  degradation  of  the  fiber  and  deterioration  of  the  fiber-matrix  interface.  This 
may  be  particularly  true  for  composites  with  glass  fibers  in  which  magnesium  and 
sodium  can  leach  out  from  the  fiber  and  possibly  contribute  to  polymer/interface 
degradation.  It  is  reported  that  graphite  fibers  containing  significant  amounts  of  sodium 
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and  potassium  as  contaminants  are  less  thermo-oxidatively  stable  than  graphite  fibers 
with  very  low  alkali  metal  contents  [18,19]. 

The  relative  oxidative  stability  of  three  polyacrylonitrile  (PAN)-based  carbon  fibers, 
namely,  T650-35,  G30-500  and  IM7,  was  recently  determined  [20]  by  comparing  the 
percent  weight  lost  during  isothermal  elevated  temperature  exposure.  As  an  example, 
Fig.  5  shows  a  comparison  of  the  weight  loss  normalized  by  the  surface  areas  of 
G30-500  carbon  fiber  and  PMR-15  neat  resin  aged  at  288°C.  Note  that  the  normalized 
fiber  weight  loss  is  almost  negligible  compared  to  that  of  neat  resin  over  2000  hrs  of 
aging.  Thus,  even  though  the  carbon  fiber  could  lose  a  significant  weight  fraction  with 
aging  at  elevated  temperature  (as  much  as  10%  for  G30-500  fibers  after  2000  hrs  of 
aging  at  288°C),  the  weight  loss  normalized  by  the  fiber  surface  area  is  negligible 
because  of  the  large  surface  area  of  bare  fibers.  Futher,  note  that  aging  of  the  bare  fibers 
may  not  necessarily  be  representative  of  the  behavior  of  the  in-situ  fibers  embedded  in 
the  matrix  because  the  exposed  surface  area  of  the  fibers  in  the  composite  is  only  a  very 
small  percentage  of  the  total  surface  area  of  the  fibers. 

While  it  has  been  established  that  the  fibers  account  for  only  a  minimal  amount  of 
the  total  weight  loss  for  HTPMCs,  nevertheless,  the  role  of  fibers  in  the  composite 
oxidation  process,  in  particular  their  role  in  transporting  oxygen  into  the  composite,  is 
not  fully  understood.  Extruded  PAN  carbon  fibers  exhibit  a  preferred  orientation  in 
which  the  graphitic  layers  may  extend  for  thousands  of  angstroms  and  extend  straight 
for  hundreds  of  angstroms  parallel  to  the  fiber  axis  providing  high-strength  and  high¬ 
stiffness  along  the  fiber  axis  direction.  It  is  very  likely  that  the  diffusivity  of  the  fibers 
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will  also  be  anisotropic  in  nature.  Unfortunately,  experimental  measures  of  carbon  fiber 
diffusivity  have  not  been  obtained.  In  addition,  while  carbon  fibers  may  be  more 
resistant  to  oxidation  and  have  longer  relaxation  times  than  the  polymer  matrix,  the 
mechanical  performance  of  the  fiber-matrix  interface  at  high  temperatures  may  be 
critical  to  the  evolution  of  damage  and  composite  failure.  The  properties  of  the  fiber- 
matrix  interface  or  interphase  region  are  a  function  of  the  interaction  of  the  fiber  and 
matrix  during  the  manufacturing  process,  the  chemical  compatibility  of  fiber,  and  the 
sizing  (if  applicable)  on  the  fiber.  The  role  of  the  fiber-matrix  interface  in  the  composite 
oxidation  process  will  be  discussed  in  subsequent  sections. 

2. 3  Composite  oxidation 

The  anisotropic  oxidative  response  of  PMCs  was  first  documented  by  Nelson  [21] 
when  he  observed  that  the  oxidation  process  is  sensitive  to  the  surface  area  for  the 
different  test  specimen  geometries  that  he  investigated.  He  found  that  the  dominant 
degradation  mechanism  for  the  graphite/polyimides  is  oxidation  of  the  matrix  at  the 
laminate  edges.  Additionally,  the  materials  degraded  preferentially  at  the  specimen 
surface  perpendicular  to  the  fiber  (axial  surface),  and  the  rate  of  oxidation  is  hastened 
by  microcracks  opening  on  the  axial  surface  increasing  the  surface  area  for  oxidation. 
The  anisotropic  nature  of  oxidation  in  HTPMCs  has  also  been  observed  by  numerous 
other  investigators  including  the  works  of  Nam  and  Seferis  [22]  and  Skontorp,  et  al.  [6]. 
Although  it  is  well  documented  that  unidirectional  composites  preferentially  oxidize  in 
the  fiber  direction  [21-23],  details  of  the  rate  of  oxidation  propagation  in  the  orthogonal 
directions  of  unidirectional  composites  is  not  well  documented. 
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Dark- field  microscopy  was  utilized  [12]  to  monitor  the  oxidation  propagation  rates 
in  both  the  axial  direction  (along  the  fiber)  and  the  transverse  direction  (transverse  to  the 
fibers)  of  unidirectional  G30-500/PMR-15  composites  aged  in  air  at  288°C.  Fig.  6 
shows  (a)  the  original  stitched  micrograph,  and  (b)  an  enhanced  micrograph  that  clearly 
distinguishes  the  white  oxidized  material  from  the  black  unoxidized  material.  The 
oxidation  layer  appears  as  a  frame  around  the  composite  specimen  just  as  seen  in  aged 
neat  resin  PMR-15  samples  in  Fig.  2.  The  specimen  of  Fig.  6  aged  for  197-hours  is 
shown  to  have  only  minimal  oxidation  transverse  to  the  fibers,  but  has  moderate 
oxidation  development  in  the  axial  direction.  The  method  of  enhancing  the  micrograph 
consists  of  constructing  (using  Adobe  Photoshop  7.0)  a  complete  image  of  the  entire 
composite  by  stitching  together  individual  micrographs  using  standard  light  microscopy 
in  the  grayscale  mode.  Once  the  image  is  constructed,  the  apparent  light  oxidized  region 
is  best  fit  in  the  lab  mode  to  a  pure  white  specified  as  having  a  lightness  value  of  100, 
while  the  remaining  unoxidized  regions  of  the  image  are  given  a  lightness  value  of  zero. 
Thus,  this  image  processing  creates  exactly  two  distinct  grayscale  colors:  black  and 
white. 

Fig.  7  shows  enhanced  micrographs  of  G30-500/PMR-15  unidirectional  composites 
after  407,  1200  and  2092  hours  of  aging  at  288°C  that  clearly  show  that  oxidation 
substantially  increases  in  both  the  axial  and  transverse  directions  with  aging  time.  Next, 
the  overall  level  of  oxidation  is  quantified  by  measuring  the  percentage  of  the  specimen 
cross-sectional  surface  area  that  is  oxidized.  Using  constructed  enhanced  micrograph 
images,  a  histogram  of  the  images  are  used  to  determine  the  ratio  of  white  to  black 
pixels,  and  the  ratios  are  used  to  quantify  the  amount  of  surface  area  of  the  composite 


cross  section  that  is  oxidized.  Quantification  of  the  oxidation  area  using  this  process  is  a 
two-dimensional  representation  and  does  not  account  for  stochastic  variations  of  the 
oxidized  cross  section  through  the  thickness  of  the  specimen.  Fig.  8  shows  the  average 
of  the  extent  of  the  axial  oxidation  along  the  fiber  direction  as  a  function  of  aging  time 
at  288°C.  The  figure  also  includes  the  average  of  the  extent  of  the  transverse  oxidation 
from  the  enhanced  micrographs.  Fig.  8  clearly  shows  the  dominance  of  the  axial 
oxidation  degradation  as  compared  to  the  transverse  degradation  in  composites  which 
can  be  attributed  to  the  higher  effective  diffusivity  on  the  axial  surfaces.  Fig.  9 
compares  the  oxidation  growth  in  the  transverse  direction  of  the  composite  to  the 
oxidation  growth  in  neat  resin  PMR-15.  It  appears  that  the  presence  of  fibers  may 
initially  retard  growth  of  oxidation  in  the  transverse  directions.  Figs.  8  and  9  clearly 
show  the  significant  differences  between  the  oxidative  behavior  of  the  constituents  as 
compared  to  the  composite  owing  to  the  fiber-matrix  microstructure  that  may  cause 
increases  in  diffusion  and  reaction  rates. 

Numerous  researchers  have  attributed  the  anisotropic  diffusivity  of  unidirectional 
laminates  to  the  preferential  oxidation  along  the  fiber-matrix  interface  or  along  the 
interphase  region1.  However,  the  interphase  region  is  a  very  small  volume  percentage  of 
the  composite  and  a  very  small  area  fraction  of  the  surface  area.  It  is  questionable  that 
the  diffusion  of  oxygen  through  and  subsequent  oxidation  of  the  interphase  region  by 
itself  can  have  such  a  significant  effect  on  the  composite  oxidation  process. 


1  Although  fiber-matrix  interface  and  fiber-matrix  interphase  are  often  used  interchangeably,  we  formally 
define  the  fiber-matrix  interface  as  the  two-dimensional  surface  defined  by  the  common  fiber  and  matrix 
surfaces.  We  define  the  fiber-matrix  interphase  as  the  three-dimensional  matrix  region  directly  around  the 
fiber  that  may  have  properties  distinct  from  the  bulk  matrix  properties. 
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More  likely,  the  rapid  oxidation  rate  along  the  fiber  length  indicates  that  in  addition 
to  diffusion  of  oxygen  from  the  specimen  surface  to  the  oxidation  front,  the  interphase 
region  has  a  supplemental  oxygen  path  or  source.  Two  possible  mechanisms  or 
scenarios  for  the  supplemental  oxygen  source  to  the  interphase  can  be  considered. 
Firstly,  if  oxidation  in  the  interphase  region  leads  to  early  fiber-matrix  interface  debond, 
the  debond  provides  a  pathway  for  oxygen  to  penetrate  deeper  into  the  composite. 
Knowing  that  there  are  residual  curing  stresses  at  the  fiber-matrix  interface,  it  is  likely 
that  such  a  scenario  would  entail  the  interface  debond  propagating  along  with  the 
oxidation  front.  The  second  scenario  is  that  the  oxygen  diffuses  into  the  composite 
along  the  fiber  and  to  the  fiber-matrix  interface  at  a  rate  much  greater  than  through  the 
neat  resin  or  interphase.  Based  on  the  anisotropic  oxidation  behavior  of  unidirectional 
composites,  the  axial  diffusivity  of  the  fibers  would  have  to  be  much  greater  than  the 
diffusivity  of  the  resin.  Unfortunately,  quantitative  measures  of  the  diffusivity  of  carbon 
fibers  are  lacking.  Although  one  might  anticipate  that  the  bulk  diffusivity  of  carbon  is 
representative  of  that  of  the  PAN  fibers,  this  would  only  be  the  case  if  the  morphology 
of  the  fibers  and  bulk  carbon  match  which  in  general  is  not  the  case. 

3.  Polymer  Oxidation  Modeling 

The  thermo-oxidation  model  for  the  polymeric  matrix  used  in  this  study  has  been 
described  in  detail  [4,5].  Thermo-oxidative  aging  was  simulated  with  a  diffusion/reaction 
model  in  which  temperature,  oxygen  concentration  and  weight  loss  effects  were 
considered.  Fig.  10  describes  the  three-region  model  in  which  the  surface  oxidative  layer  is 
separated  from  the  unoxidized  polymer  with  an  active  reaction  zone.  As  oxygen  diffuses 
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through  the  polymer,  it  is  consumed  by  the  oxidation  reaction.  Once  a  region  is  fully 
oxidized,  oxygen  can  diffuse  through  it  and  oxidation  begins  in  the  adjoining  region.  This 
creates  three  distinct  zones.  Closest  to  the  oxygen  source  is  a  region  of  completely 
oxidized  polymer.  This  is  followed  by  the  transition  region,  where  oxidation  is  actively 
occurring,  and  finally  by  the  unoxidized  interior  of  the  material.  Although  this  three-phase 
description  is  analogous  to  the  two-phase  models  (or  skin/core  models)  described  by  Nam 
and  Seferis  [22],  the  description  of  three  phases  enables  extension  to  situations  with  high 
conversion  ratios  in  which  the  constancy  of  the  substrate  concentration  is  no  longer  valid 
and  the  substrate  depletion  leads  to  an  auto-retardation  of  the  oxidation  kinetics.  A 
practical  consequence  is  a  decrease  of  the  oxidation  rate  in  the  superficial  (region  I)  layers 
of  the  specimen,  where  the  high  conversion  ratios  are  firstly  reached,  which  leads  to  a 
moving  of  the  oxidation  front  towards  the  specimen  core.  Such  an  extension  was  proposed 
by  Colin,  et  al.  [24,25]  and  successfully  applied  to  the  thermal  oxidation  of  amine  cross- 
linked  epoxy  and  bismaleimide  resins. 

The  molar  concentration  of  oxygen  at  any  material  point  of  a  constituent 
material,  C(x,t ) ,  is  determined  by  the  diffusion  reaction  equation  (Eq.  1). 

p\<t>,T)\2C(x,t)-R\C)  (1) 

sf5  Sfc 

where  D  is  the  diffusivity  and  R  (C)  represents  the  oxidation  reaction  rate  function  in  the 
material.  The  exposed  boundaries  are  subject  to  the  oxygen  sorption  as  governed  by  the 
Henry’s  Law.  The  exposed  boundaries  are  therefore  assigned  a  molar  concentration  value 
equal  to  C  and  oxygen  flow  across  symmetry  boundaries  is  set  to  zero  (Eq.  2).  A 


dC(x,t ) 
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negligible  initial  concentration  (s  ~le-7  mol/m3)  is  assumed  through  out  the  domain  to 
avoid  numerical  difficulties. 

C  =  Cs  on  the  exposed  boundaries 
VC.  n  =  0  on  symmetry  boundaries  (2) 

C(x,0)  =  e 

We  introduce  an  oxidation  state  parameter,  (f>,  which  tracks  the  progress  of  the 
oxidation  reaction.  <fi  varies  from  (j)  =  1  for  unoxidized  state  to  (j)  =  </>ox  for  completely 
oxidized  state.  The  determination  of  <f>ox  value  is  explained  in  earlier  papers  [4,5].  The 
reaction  rate  function  R*(C)  is  influenced  by  the  temperature  T,  available  oxygen 
concentration  C,  and  is  also  a  function  of  the  polymer  substrate  oxidation  state  variable 
(f>  as  shown  in  Eq.  3 

R*(C)  =  g(</>)R(C,T)  (3) 

where  g(<f>)  describes  the  reaction  rate  dependence  on  the  oxidation  state  variable,  and 
R{C,T )  defines  the  reaction  rate  dependence  on  the  concentration  and  temperature.  The 
reaction  kinetics  are  modeled  following  Colin,  et  al.  [24,25]  and  Abdeljaoued  [11]  in 
which  a  constant  saturation  reaction  rate  R0  is  the  maximum  reaction  rate  when  reaction 
is  not  oxygen  starved  (Eq.  4).  The  saturation  reaction  rate  is  temperature  dependent  and 
an  Arrhenius  relationship  is  typically  assumed  for  the  temperature  dependence, 
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where  the  parameter  f>  nondimensionalizes  the  concentration  field,  Rq  is  the  rate  constant, 

and  Ra  is  the  activation  parameter.  Eq.  5  shows  the  diffusivity  of  an  oxidizing  constituent. 
The  diffusivity  will  depend  upon  the  oxidation  state  and  temperature.  Therefore,  two 
values  of  diffusivity  (assuming  isotropic  diffusivity,  scalar  diffusivity  value  is  represented 
instead  of  the  tensor  values  representing  anisotropy)  for  the  oxidized  and  unoxidized  states 
of  the  constituent  are  denoted  by  ox  and  un  superscripts,  respectively. 


DU,  T )  =  Dl,n  (7)f -  +  Dox  (7)1  1  ^ 

IWoJ  V  ~<t>ox 

1 7un 
~ha 

Dun  (T)  =  D"n  e  RT  (5) 

_FOX 

Dox{T)  =  D°x  e  RT 

The  oxidation  state  parameter  can  be  directly  correlated  with  the  oxidative  weight  loss 
(weight  loss  due  to  oxidation  alone,  typically  obtained  by  subtracting  the  weight  loss  in  an 
inert  -  argon  -  environment  from  that  in  an  oxidative  -  air/oxygen  -  environment).  A  time 
dependent  proportionality  relationship  is  assumed  between  the  reaction  rate  and  the  rate  of 
oxidation  as  shown  in  Eq.  (6).  The  reaction  rate  is  assumed  to  have  a  linear  dependence 
during  the  first  tswitch  hours  of  aging,  reducing  from  amax  to  a  value  of  arain  within  tswitch 
hours.  The  reaction  rate  after  tswitch  hours  is  assumed  to  be  constant  and  equal  to  lower  amin. 
Again,  we  refer  you  to  references  [4,5]  for  further  details  on  the  rationale  of  this  selection 
and  correlations  with  the  experimental  results. 
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d(j) 

dt 


-a(t)R(C,T) 


a(t)= 


„  max  „  max  „  mm 

a  -\a  -a 


t 


\  t  switch 

min  .  .  . 

CC  .  I  >  l switch 


■  t  <  t switch 


(6) 


The  diffusion-reaction  system  of  equations  (1-6)  can  be  solved  using  numerical 
solutions  to  differential  algebraic  equations.  The  capability  to  obtain  direct 
measurements  of  all  of  the  properties  of  the  PMR-15  for  each  of  the  three  material 
regions  that  are  needed  to  populate  the  parameters  for  the  model  is  lacking.  Therefore, 
many  of  the  parameters  were  obtained  [4,5]  through  indirect  measurement  or  correlation 
with  model  predictions  of  the  oxidation  layer  growth  obtained  at  288°C. 

Fig.  1 1  shows  predictions  of  the  oxidation  layer  thicknesses  for  288°C  (550°F), 
316°C  (600°F),  and  343°C  (650°F)  neat  resin  PMR-15  specimens  based  on  288°C 
(550°F)  properties  shown  in  Tables  1-3.  The  figure  shows  that  the  model  accurately 
predicts  the  316°C  (600°F)  and  343°C  (650°F)  oxidation  layer  thicknesses,  and  the 
temperature  dependence  of  the  reaction  rate  and  diffusivity  are  accurately  represented. 
Further,  as  mentioned  earlier,  acceleration  of  oxidation  growth  can  be  achieved  by 
increasing  the  partial  pressure  of  oxygen  within  the  aging  chamber.  Since  the  partial 
pressure  of  a  gas  is  directly  proportional  to  the  total  pressure  and  in  general  the 
solubility  of  gases  increase  with  increasing  pressure,  the  sorption  is  expected  to  increase 
with  increasing  pressure.  Fig.  12  shows  the  results  of  a  parametric  evaluation  of  how 
changes  in  the  sorption  affect  the  predicted  oxidation  layer  thickness.  Oxidation  layer 
thickness  for  neat  resin  PMR-15  specimens  aged  at  ambient  and  elevated  pressures  at 
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288°C  are  also  shown  in  the  figure.  The  predicted  oxidation  layer  thickness  increases 
with  increasing  oxygen  concentration  on  the  boundary  and  for  the  sorption  value 
Cs  =  1.25  ,  there  is  a  good  correlation  with  the  pressurized  aging  data  for  aging  times  up 
to  200  hours.  Fig.  13  shows  the  predicted  oxidation  layer  thickness  for  predictions  up  to 
800  hours  of  aging  time.  The  results  show  good  agreement  with  the  experimental  data. 

4.0  Composite  Oxidation  Modeling 

4.1  Finite  Element  Framework  for  Modeling  Thermo-Oxidation  in  Composites 

Both  two-dimensional  and  three-dimensional  unit  cell  models  capable  of  discrete 
crack  analysis  coupled  with  thermo-oxidative  aged  analysis  can  be  used  to  model  and 
identify  the  origins  of  degradation  and  mechanisms  of  failure  initiation  representative  of 
that  in  unidirectional  composites.  Traditionally  the  fibers  are  assumed  to  be  chemically 
static  for  short  aging  times  (this  assumption  may  not  be  valid  for  extended  aging  times) 
and  residual  stresses  must  be  taken  into  account  due  to  the  strong  coupling  between 
stress,  diffusion,  and  thus  aging.  Of  primary  importance  is  the  characterization  of  the 
interface/interphase  behavior  that  may  have  different  reaction/diffusion  kinetics  than  the 
polymer  matrix  phase  of  the  composite.  The  fiber-matrix  interphase  region  developed 
during  the  composite  cure  cycle  is  likely  dependent  on  the  fiber  surface  (unsized  or 
sized)  and  the  composition  of  the  fiber  sizing,  if  it  is  present. 

The  diffusion-reaction  oxidation  model  is  implemented  into  a  three-dimensional 
Galerkin  Finite  Element  Analysis  (GFEA).  The  GFEA  method  is  more  appropriate  for 
studying  structural  scale  problems  than  the  one-  and  two-dimensional  implementations; 
however,  the  computational  requirement  of  the  three-dimensional  analysis  precludes  its 
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use  in  parametric  analyses  required  for  correlation  with  experimental  results  and 
parameter  determination.  The  GFEA  requires  mesh  sizes  in  the  1  -pm  scale  and  time 
increments  in  1 -minute  steps.  A  200-hour  oxidation  simulation  with  100-pm  oxidation 
zone  size  typically  requires  problem  sizes  in  the  order  of  500,000  degrees  of  freedom 
(DOF)  and  12,000  time  steps.  The  domain,  Q ,  is  discretized  by  a  union  of  finite 

elements  as  Q  =  [J N‘_  ,  where  =  0 ,  if  j ,  and  Ne  is  the  total  number  of 

elements  in  Q .  The  approximate  element  solution  for  the  concentration  C (x,y,z,t)  in  the 
diffusion  model  is  defined  as 

C  (x,  y,  z,  t)  =  [®  (x,  y,  z)  f  [c{t)\  (7) 

where  O (x,y,z)  is  the  spatial  interpolation  function  and  C(t)  is  the  time  varying  nodal 
concentration  values.  Using  the  Galerkin  method,  i.e.  using  a  weight  function  same  as 
the  interpolation  O , 

w  =  Q>(x,y,z)  (8) 

the  residual  for  the  diffusion-reaction  equation  (Eq.  7)  must  vanish  as  given  in  Eq.  9 

j)[^0®r  ^llyv  +  ^R\C))dV 

rf  _ .  f  50  50  r  ,  ,  f 50  50 r  .  f 50  50 T  .  vYl  . T/  A 

+  J  Al  ~  ~  A2  ~  ~  ^33  ~  dV  —  0  (  ) 

y  V  vx  dx  y  y  oy  oy  y  y  oz  oz  yy 

where  orthotropic  diffusivity  is  assumed.  The  superscript  *  on  the  diffusivity  tensor,  D, 

sjc 

and  reaction  rate,  R  (C)  is  used  to  indicate  that  the  diffusivity  and  reaction  rate 
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functions  are  dependent  on  the  temperature  T  and  oxidation  state  variable,  </>  .  After 


some  mathematical  manipulations,  Eq.  9  can  be  written  in  matrix  form  as: 


[c]5^)  +  [b]c(0+{r}  =  o 

at 


(10) 


where,  the  matrices,  [B],  [C]  ,  and  vector  {r}  are  given  by, 
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Eq.  10  is  discretized  by  backward  Euler  method  in  time  domain,  where  At  =  tn+1  - 1" . 
Therefore,  Eq.  10  can  be  denoted  as  a  system  of  linear  algebraic  equations 


[C} 


VL 

At 


[b]|c"+1  +{r}  =  o 


(14) 


In  summary,  the  diffusion-reaction  equations  in  their  discretized  form  can  be  written 
as: 

<R  =  [d]{c"+1}-{e}  =  0  (15) 

where, 


[d]=[c]+a<[b] 

(16) 

[e]  =  [c]{c"  |-  {r} 

(17) 
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Eq.  15  is  solved  using  a  direct  solution  or  by  minimizing  the  residual  ( )  using  an 
iterative  solution  such  as  Jacobi  Conjugate  Gradient  (JCG)  iterative  method  [26]  for 
problems  with  a  large  number  of  degrees  of  freedom.  The  time  step,  At,  used  for  this 
parametric  study  is  1  minute  and  the  residual  is  minimized  to  a  root  mean  square  (L2 
norm)  value  of  lxlO"8. 

A  simple  user  interface  for  three-dimensional  thermo-oxidative  modeling  labeled 
TOM-3D  is  established  to  set  up  the  problem  domains,  material  properties  and  the 
oxidation  rates  for  each  material  as  shown  in  Fig.  14.  TOM-3D  first  calculates  the 
concentration  profiles  and  then  updates  the  oxidation  state  (4>)  for  each  element.  The 
oxidation  state  is  taken  to  be  constant  over  each  time  step  and  used  for  calculation  of  the 
(^-dependent  diffusivity  and  reaction  rates. 

Typical  meshes  used  for  composite  unit  cells  are  shown  in  Fig.  15.  Fig.  15(a)  shows 
a  single  fiber  representative  volume  element  (RVE)  with  the  fiber  end  exposed  to 
oxidation  to  simulate  axial  oxidation  growth.  It  is  assumed  that  the  fibers  are  in  a  square 
arrangement,  with  each  fiber  having  a  diameter  of  10  pm  and  a  spacing  of  12.532  pm 
(center  to  center)  between  adjacent  fibers.  The  RVE  is  taken  to  have  a  nominal  fiber 
volume  fraction  of  50%,  but  can  be  varied  easily  by  changing  the  FEM  mesh.  The 
length  of  the  model  is  scaled  appropriately  to  the  expected  oxidation  size.  The  shortest 
model  had  a  length  of  125  pm  and  the  longest  1000  pm  in  this  parametric  study. 

Though  not  presented  in  this  paper,  we  have  further  observed  that  the  simulations  for 
oxidation  growth  in  axial  direction  using  a  single  fiber  model  closely  match  those 
obtained  using  multiple  fibers  with  both  the  fiber  ends  and  axial  surface  parallel  to  the 
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fibers  exposed  to  oxidation.  Further,  edge  effects  for  fibers  located  close  to  the  exposed 
axial  surface  dissipate  within  a  short  time  period  (~  40  hours).  Since  the  transverse 
oxidation  for  the  composite  can  span  multiple  fibers,  a  ten  fiber  125  pm  x  12.5  pm  x 
125  pm  (shown  in  Fig.  15(b))  RVE  and  a  20  fiber  250  pm  x  12.5  pm  x  125  pm  RVE 
are  created  for  the  transverse  simulations.  The  unit-cell  size  along  the  length  of  the  fiber 
(axial  direction)  is  scaled  to  accommodate  the  estimated  oxidation  layer  and  reaction 
zone  sizes.  Note  that  for  this  unit  cell,  both  axial  and  transverse  boundaries  are  exposed 
to  oxidation.  Similar  to  axial  oxidation  growth,  edge  effects  in  transverse  simulation 
become  insignificant  in  a  relatively  short  time  period.  An  interphase  region  of  arbitrary 
thickness  has  also  been  modeled  to  create  a  third  phase  between  the  fiber  and  matrix 
domains.  The  size  and  material  properties  of  the  interphase  region  can  be  parametrically 
varied.  The  10  fiber  mesh  used  has  352,012  elements  with  481,761  nodal  degrees  of 
freedom  (DOF)  for  concentration  and  same  number  for  the  oxidation  state  for  a  total  of 
963,522  DOFs.  The  20  fiber  RVE  has  447,134  DOFs  and  154,486  elements.  The  20 
fiber  RVE  is  coarser  in  comparison  to  10  fber  RVE  in  order  to  minimize  the 
computational  time.  Also,  the  reaction  zone  size  for  the  20  fiber  model  is  expected  to  be 
larger  because  of  the  selected  model  parameters,  thus,  mitigating  the  need  for  a  finer 
mesh,  as  required  for  the  1 0  fiber  simulation.  In  the  analyses  presented  in  this  paper,  it 
is  assumed  that  the  fibers  do  not  react  with  oxygen  and  are  purely  diffusive,  the 
interphase  diffusivity  is  parametrically  varied  (as  it  is  problematic  to  directly  measure 
interphase  diffusivity),  and  oxidation  only  occurs  in  the  resin.  We  also  set  the  sorption 
value  of  the  fiber  and  interphase  regions  to  be  the  same  as  that  of  the  matrix  (Cs  =  0.79 
mol/m3).  Each  of  the  constituents  is  modeled  explicitly  as  a  separate  material  and  the 
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interface  surfaces  at  the  fiber/interphase  and  interphase/resin  boundaries  are  assumed  to 
be  perfectly  diffusive.  Resin  properties  corresponding  to  that  of  PMR-15  as  given  in 
Tables  1-3  are  used. 

4.2  Composite  Oxidation  Studies  with  Representative  Volume  Elements  (RVE’s) 

In  this  section,  the  diffusivity  of  the  fiber  is  parametrically  varied  to  evaluate  its 
influence  on  oxidative  growth.  Let  Df,  1),  and  Dm  denote  the  constituent  level  isotropic 
diffusivity  values  of  the  fiber,  interphase  region  and  matrix,  respectively.  A  three- 
dimensional  RVE  is  used  with  PMR-15  matrix  and  an  interphase  region  is  initially  not 
included.  The  oxidation  simulation  is  mn  for  up  to  200  hours  because  steady  state 
oxidation  layer  growth  can  be  expected  at  that  time.  The  oxidation  state  in  the  model  is 
plotted  at  various  time  intervals  (25,  50,  100  and  200  hours).  Fig.  16  shows  the 
oxidation  growth  in  both  axial  and  transverse  directions  for  various  fiber  diffusivities 
and  at  several  hours  of  exposure.  Let  the  oxidation  thickness  in  the  axial  and  transverse 
direction  be  denoted  by  taxiai  and  tlrans,  respectively.  Similarly,  we  will  denote  the 
thickness  of  the  reaction  zone  in  the  axial  and  transverse  direction  by  laxiai  and  ltrans> 
respectively.  The  oxidation  layer  sizes  (tax;ai  and  ttrans)  and  the  reaction  zone  size  (laxjai 
and  brans)  are  measured  using  the  oxidation  state  (4»)  values.  The  oxidized  layer  has  been 
considered  to  have  4>  <  0.181  (99.5%  oxidation)  and  the  unoxidized  layer  with  4>  > 
0.998  (corresponding  to  0.2  %  oxidation  in  the  material).  These  limits  are  primarily 
chosen  to  ensure  consistent  measurement  of  the  zone  sizes  throughout  the  parametric 
study. 
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Fig.  16(a)  shows  the  case  in  which  fiber  has  a  very  low  diffiisivity,  equivalent  to 
one  hundredth  that  of  the  oxidized  resin,  Fig.  16  (b)  corresponds  to  that  with  a  fiber 
diffusivity  of  one  tenth  that  of  the  resin  and  Fig.  16  (c)  shows  the  oxidation  layer  size 
when  the  fiber  diffusivity  is  about  ten  times  that  of  the  resin.  For  these  simulations,  only 
the  diffusive  behavior  of  the  fibers  is  varied  whereas  both  the  diffusion  and  reaction 
behavior  of  the  resin  is  considered.  When  the  oxidation  state  of  the  resin  reaches  full 
oxidation  (<j)=(j)ox=0.18),  the  oxidized  resin  becomes  non-reactive  and  diffusive  only. 
Prior  analysis  shows  that  the  diffusivity  of  the  oxidized  matrix  is  higher  than  that  of  the 
unoxidized  matrix. 

When  the  fiber  is  less  diffusive  than  the  matrix  as  in  cases  16  (a)  and  16  (b),  the 
oxidation  layer  size  in  both  axial  and  transverse  directions  is  not  significantly  influenced 
by  the  fiber  diffusivity.  Thus,  for  low  fiber  diffusivity  values,  oxygen  diffusion  through 
the  fiber  becomes  insignificant  and  oxidation  is  driven  by  matrix  diffusivity  alone, 
ffowever,  when  the  diffusivity  of  the  fiber  is  higher  than  that  of  the  matrix  as  in  Case  16 
(c),  not  only  is  the  oxidation  layer  size  larger,  the  active  reaction  zone  size  is  also  larger 
indicating  the  additional  supply  of  oxygen  to  the  active  oxidation  region  of  the 
composite.  The  effective  diffusivity  of  the  oxidized  region  of  the  composite  controls 
the  growth  rate  and  extent  of  the  oxidation  layer. 

It  is  generally  believed  that  the  oxidation  process  in  PMR-15  resin  is  diffusion 
limited  versus  reaction  rate  limited.  Although  the  constituent  level  diffusivities  of  the 
composite  are  assumed  to  be  isotropic,  the  geometry  of  the  unidirectional  composite 
results  in  effective  diffusivities  that  are  transversely  isotropic.  Thus,  in  order  to  establish 
the  relationship  between  the  homogenized  diffusivities  of  the  oxidized  composite  and 
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the  observed  oxidation  layer  size,  we  calculate  the  ratio  of  the  oxidation  layer  size  in  the 
composite  in  the  axial  direction  (taxiai)  to  that  in  the  transverse  direction  (ttrans)  and  as  a 
first  approximation  compare  it  with  the  corresponding  ratio  of  homogenized 
diffusivities  of  the  composite.  The  composite  effective  diffusivities  were  recently 
calculated  [27]  using  homogenization  theories  for  various  ratios  of  Dt/Dm  and  are  listed 
in  Table  4.  Here  Di  and  D2  denote  the  effective  axial  and  transverse  diffusivities, 
respectively,  of  the  unidirectional  composite.  Table  4  also  shows  a  comparison  of  the 
anisotropy  ratio  of  the  diffusivities  (D1/D2)  with  the  corresponding  oxidation  layer  ratio 
(taxiai/hrans)  and  the  reaction  zone  ratio  (laxiai/ltrans)  at  200  hours  simulations  for  some 
selected  fiber  and  matrix  diffusivity  values.  For  values  of  Dt/Dm  <  10,  the  oxidation 
growth  is  diffusion  limited  and  the  diffusivity  factor  obtained  from  the  homogenization 
method  closely  matches  the  observed  anisotropy  in  the  oxidation  layer  sizes.  However 
for  larger  ratios  of  Dt/Dm  (i.e.  ,>  10),  the  oxidation  process  is  observed  to  become 
reaction  rate  limited  as  the  reaction  zone  ratio  also  increases  with  increasing  Dt/Dm 
ratio.  The  oxidation  layer  ratio  taxiai/ttrans  therefore  is  calculated  to  be  less  than  the 
corresponding  D1/D2  diffusivity  factor  since  oxidation  is  not  occurring  as  rapidly  as  the 
rate  at  which  oxygen  is  being  supplied.  Hence  a  larger  portion  of  the  surface  layer  now 
has  regions  in  which  complete  oxidation  has  not  occurred,  leading  to  an  increase  in  the 
size  of  the  “active”  reaction  zone. 

In  earlier  discussion,  we  had  suggested  that  the  diffusivity  of  carbon  fibers  will 
likely  be  anisotropic  in  nature  since  the  extruded  PAN  carbon  fibers  exhibit  a  preferred 
orientation.  To  examine  this,  we  simulated  the  oxidation  growth  with  several  assumed 
values  of  fiber  axial  and  transverse  diffusivity  values  denoted  by  Df, axial  and  Dfjtrans 
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respectively,  and  compared  the  solution  with  the  corresponding  isotropic  fiber  case. 
Table  5  lists  the  results  of  our  simulations.  As  an  example,  an  isotropic  fiber  with  Df  = 
0.1  Dmox  produces  a  taxiai/ttrans  value  of  1.31  (taxiai=  57.5  pm  ,  ttrans=  43.9  pm)  and  with 
Df  =  10  Dmox  produces  1.56  (taxiai=  122.5  pm  ,  ttrans=  78.56  pm).  In  comparsison,  an 
anisotropic  fiber  with  Dr,axsai  =10  Draox  and  Doans  =0. 1  Dmox  produces  a  taxiai/ttrans  value 
of  4.41  after  200  hours  of  oxidation  (tax;ai=  193.75  pm  ,  t,rans  =  43.9  pm).  Thus,  the 
simulations  show  that  anisotropy  in  fiber  diffusivity  can  lead  to  larger  oxidation  layer 
ratio  in  unidirectional  composites.  Similar  trends  are  observed  for  other  diffusivity 
values  as  seen  in  Table  5.  However,  the  extent  of  anisotropy  seen  in  the  simulations  is 
not  as  much  as  that  seen  in  experimental  observations  of  oxidation  growth  in 
G30-500/PMR-15  composites  where  the  ratio  (taxiai+laxiai)/(ttrans+ltrans)  approaches  nearly 
eight. 

4.3  RVE  with  Interphase  Layers 

To  investigate  the  role  of  the  interphase  region,  particularly  if  the  interphase 
thickness  and  its  diffusivity  are  large,  a  study  is  conducted  to  understand  the  influence 
of  the  interphase  parameters  on  the  oxygen  diffusivity  and  transport  in  the  composite. 
This  is  accomplished  by  extending  the  RVE  analysis  to  include  a  finite  thickness 
interphase  layer.  The  interphase  layer  is  assumed  to  be  a  concentric  cylinder  with  that  of 
the  fiber  with  a  thickness  of  5%  that  of  the  fiber  radius.  We  assume  the  fiber  diffusivity 
to  be  isotropic  and  set  it  equal  to  0.1  Dmox,  while  the  diffusivity  of  the  interphase  is 
assumed  to  be  anisotropic,  with  axial  diffusivity  f)Laxiai  =  10  Draox  and  transverse 
diffusivity  Dorans  =  Dmox.  Fig.  17  shows  the  comparison  of  the  oxidation  layer  growth 
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within  an  RVE  which  includes  the  interphase  with  another  RVE  in  which  no  interphase 
is  present  at  various  time  intervals  (10,  20  and  30  hours).  Clearly,  the  presence  of  a 
highly  diffusive  interphase  region  contributes  to  an  increase  in  the  oxidation  layer 
growth  transverse  to  and  in  the  fiber  direction.  The  evidence  of  higher  diffusivity  in  the 
axial  direction  compared  to  that  in  transverse  direction  can  also  be  found  in  the  sizes  of 
the  active  reaction  zones.  Next,  for  results  shown  in  Table  6,  we  set  the  fiber  diffusivity 
at  0.1  Dmox  while  the  interphase  diffusivity  values  are  varied.  The  first  row  in  Table  6  is 
the  result  at  200  hours  simulation  with  no  interphase  region.  We  observe  that  the 
addition  of  a  highly  diffusive  anisotropic  interphase  results  in  substantial  increase  in 
composite  oxidation  in  the  axial  direction,  while  only  small  changes  occur  in  oxidation 
region  in  the  transverse  direction.  Thus,  similar  to  the  fiber  anisotropy,  the  simulations 
show  that  anisotropy  in  interphase  diffusivity  can  also  lead  to  larger  oxidation  layer 
ratio  in  unidirectional  composites. 

4. 4  Modeling  with  Damage  Evolution  at  the  Interphase 

Our  modeling  results  so  far  indicate  that  the  experimentally  observed  extent  of 
anisotropy  in  the  oxidation  growth  cannot  be  explained  with  rational  choices  for  fiber 
and  interphase  diffusivities  alone.  The  rapid  oxidation  rate  along  the  fiber  length 
indicates  that  a  degradation  process  is  creating  oxygen  diffusion  paths  deep  into  the 
specimen,  possibly  in  the  form  of  fiber-matrix  debonding.  In  order  to  simulate  the 
effect  of  interphase  damage  on  the  oxidation  layer  size,  an  oxidation  state  dependent 
damage  evolution  model  can  be  formulated  by  extending  the  basic  framework  of  the 
matrix  oxidation  model.  The  damage  evolution  itself  can  be  considered  to  be  similar  to 
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that  of  the  oxidation  reaction  and  consumption  model  where  in  the  interphase’s  damage 
state  can  be  tracked  by  a  damage  state  parameter  (4>d).  The  damage  process  can  be 
controlled  by  the  presence  of  oxygen  (through  the  oxygen  concentration  variable, 
c(x,y,z,t))  and  a  reaction  rate  process,  R,(C).  The  damage  process  can  therefore  be 
analogous  to  subtrate  conversion  in  the  oxidation  simulation  and  the  rate  can  be 
controlled  by  a  conversion  parameter,  a  a  which  will  be  analogous  to  the  subtrate 
conversion  rate  parameter  a.  Next,  a  diffusivity  to  damage  relation  will  have  to  be 
assumed  for  the  interphase  region.  As  the  damage  state  of  the  material  will  change  from 
undamaged  to  damaged,  the  diffusivity  of  the  interphase  can  be  altered  in  the 
simulation.  This  will  allow  the  damage  and  the  thermo-oxidation  processes  to  follow 
independent  rate  kinetics  as  we  simulate  the  effect  of  damage  on  the  thermo-oxidation. 

5.  Summary 

Thermo-oxidative  degradation  of  PMCs  is  clearly  a  subject  of  major  importance  as 
designers  push  the  high  temperature  endurance  limits  of  polymer  composites  to  improve 
the  performance  of  aerospace  systems.  Along  with  all  of  the  tremendous  challenges  of 
modeling  the  performance  and  strength  of  composites  in  benign  thermal  environments, 
are  added  the  difficult  challenges  of  modeling  the  age-dependent  properties  of 
composites  at  high  temperature.  The  strong  interaction  of  the  physical,  chemical,  and 
mechanical  aging  behavior  of  HTPMCs  is  due  to  the  significant  contrast  in  the 
properties  of  the  fiber  and  matrix  constituents  of  the  composite.  These  interactions  are 
particularly  apparent  in  the  critical  fiber-matrix  interface/interphase  where  the 
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interface/interphase  has  been  shown  to  play  a  dominant  role  in  the  thermo-oxidation 


process. 

In  this  work,  an  overview  of  the  experimentally  observed  thermo-oxidative  behavior 
of  PMR-15  neat  polymer  and  carbon  fiber/PMR-15  composites  has  been  presented 
along  with  a  theoretical  treatment  of  oxidation  development.  Literature  search  indicates 
that  although  the  degradation  of  the  in-situ  fiber  in  the  composite  has  long  been 
considered  to  be  negligible,  the  role  of  the  fiber  and  its  contribution  to  the  oxidation 
process  has  not  been  fully  explored.  Specifically,  little  is  known  about  the  diffusivity  of 
oxygen  in  the  fiber  and  its  contribution  to  the  transport  of  oxygen  to  the  interior  of  the 
composite.  A  method  for  high  resolution  oxidation  modeling  capable  of  simulating 
thermo-oxidation  for  composites  with  multiple  phases  and  anisotropic  diffusivity  has 
been  developed.  The  new  method  utilizes  a  diffusion/reaction  model  in  which 
temperature,  oxygen  concentration,  and  weight  loss  effects  are  considered.  Three 
dimensional  FEM  analyses  have  been  used  to  demonstrate  the  interactions  between 
fiber  and  resin  in  PMCs.  Parametric  studies  illustrating  the  anisotropy  in  the  oxidative 
region  growth  and  the  effect  of  fiber  and  interphase  diffusivity  on  the  oxidation  layer 
growth  rates  are  discussed.  It  is  suggested  that  alternative  pathways  for  transport  of 
oxygen  into  the  interior  of  the  composite  are  fiber-matrix  debonds  that  propagate  with 
the  oxidation  front.  This  possible  mechanism  for  accelerated  oxidation  along  the 
composite  fiber  direction  is  an  excellent  example  of  the  intrinsic  coupling  of  chemical 
oxidative  aging  and  damage.  Whereas  cracks  can  accelerate  oxidation  by  providing 
pathways  for  oxidants,  so  can  oxidation  lead  to  cracking  in  the  resin  matrix  and  at  fiber- 
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matrix  interfaces.  The  critical  nature  of  the  fiber-matrix  interphase  on  degradation  and 
failure  processes  in  composites  signifies  the  importance  of  proper  representation  of  its 
behavior  in  predictive  models.  Since  failure  initiation  is  typically  associated  with  the 
fiber-matrix  interface/interphase  region,  the  importance  of  modeling  the  coupling 
effects  of  damage  and  oxidation  cannot  be  over  emphasized. 

Simulations  have  compared  well  with  experimental  data,  but  further  refinement  is 
still  necessary.  While  it  is  only  one  mechanism  in  the  aging  process  of  PMCs, 
successful  simulation  of  thermo-oxidation  behavior  is  a  critical  part  of  understanding 
this  complex  process  and  in  providing  a  foundation  for  understanding  the  mechanisms 
of  degradation  for  long-term  environmental  exposure  and  loading. 
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Figure  Captions 

Fig.  1.  Oxidation  measurement  procedures 

Fig.  2.  Photomicrograph  of  PMR-15  resin  showing  the  formation  of  oxidized  region 
Fig.  3.  Evolution  of  oxidation  layer  and  transition  region  thickness  with  aging  time  in 
PMR-15  resin  at  288°C 

Fig.  4.  Oxidation  layer  thickness  for  ambient  and  elevated  pressure  aging  in  PMR-15 
Fig.  5.  Normalized  weight  loss  of  G30-500  fiber  and  PMR-15  resin 
Fig.  6.  (a)  Array  of  eight  stitched  photomicrographs  of  oxidized  specimen  cross 
section,  (b)  Enhanced  micrograph  of  the  specimen  cross  section 
Fig.  7.  Unidirectional  G30-500/PMR-15  composite  at  407,  1200  and  2092  hours  of 
oxidation  at  288°C 

Fig.  8.  Comparison  of  oxidation  growth  in  axial  and  transverse  directions  of  composite 
Fig.  9.  Comparison  of  oxidation  growth  in  transverse  direction  of  composite  and  in 
neat  resin  PMR-15 

Fig.  10.  Schematic  of  the  three-zones  in  thermo-oxidation. 

Fig.  11.  Prediction  of  oxidation  thickness  at  316°C  and  343°C  based  on  288°C  model 
parameters 

Fig.  12.  Influence  of  sorption  on  the  oxidation  layer  thickness 

Fig.  13.  Prediction  of  oxidation  thickness  for  pressurized  aging  of  PMR-15 

Fig.  14.  User  input  interface  for  TOM-3D 

Fig.  15.  Representative  volume  elements  used  for  the  analyses  (a)  Single  Fiber  RVE 

(12.5  pm  x  12.5  pm  x  125  pm)  with  one  surface  exposed  to  oxidation;  (b)  Ten 
Fiber  RVE  (125  pm  x  12.5  pm  x  125  pm)  with  both  axial  and  transverse 
directions  exposed  to  oxidation. 

Fig.  16:  Effect  of  fiber  diffusivity  on  composite  oxidation.  Black  regions  are  completely 
oxidized  and  grey  regions  are  unoxidized.  The  transition  region  is  contoured 
from  black  to  grey  based  on  the  oxidation  state.  All  dimensions  shown  are  in 
pm.  (a)  Df  =  0.01  Dmox,  (b)  Df=  0.1  Dmox,  (c)  Df=  10  Draox 
Fig.  17.  Effect  of  a  highly  diffusive  fiber-matrix  interphase  region.  The  anisotropy  in 
the  oxidation  layer  can  be  clearly  seen. 
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Table  1 


Diffusivity  parameters  for  PMR-15  resin 


Unoxidized,  (/>  =  1 

Oxidized,  (j)  =  c/)ox 

Do 

8.90x1 011  m2/sec 

Ea 

19,700  J/mol 

R 

8.31447  J/(mol  °K) 

Table  2 

Reaction  rate  parameters  for  PMR-15  resin 


Parameter 

Value 

Ro 

536.94  mol/(mJ.min) 

Ra 

23379.63  J/mol 

P 

Table  3 


Weight  loss  model  parameters  for  PMR-15  resin 


Parameter 

Value 

a 

0.01-(0.01-.0033)*t/40  mJ/mol,  t<40hrs 
0.0033  m3/mol,  t  >  40  hrs 

tyox 

0.18 
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Table  4 


Anisotropy  in  composite  diffusivity  and  its  correlation  with  oxidation  layer  size 


Diffusivity /Dmox 

Diffusivity  Factor 

Predicted 
Oxidation 
Layer  Size 
Ratio 

(taxial/ttrans) 

Predicted 
Reaction 
Zone  Size 
Ratio 

Oaxia|/f  rails) 

Matrix 

Fiber 

D,/Dmox 

D2/Dmox 

Di/D2 

1 

0.1 

0.551 

0.415 

1.33 

1.31 

1.15 

2 

1.499 

1.4 

1.07 

1.20 

0.83 

4 

2.34 

1.653 

1.42 

1.43 

0.95 

10 

5.49 

2.41 

2.28 

1.56 

2.00 

20 

10.48 

2.69 

3.89 

2.71 

2.68 

46 

23.45 

2.89 

8.11 

3.25 

3.5 

60 

30.52 

2.91 

10.49 

3.45 

3.93 

100 

50.37 

2.99 

16.84 

4.41 

5.39 

200 

100.5 

3.015 

33.33 

7.40 

10.53 
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Table  5 


Influence  of  fiber  anisotropy  on  composite  oxidation 


Df,  axial / 

D  ox 

Df,  trans/ 

D  ox 

taxial, 

Jim 

laxial, 

jam 

ttrans, 

Jim 

ltrans, 

Jim 

taxial/ 

ttrans 

(taxial  laxial)/ 
(ttrans-*-  ltrans) 

0.1 

0.1 

57.5 

23.75 

43.9 

20.6 

1.31 

1.15 

1.26 

10 

10 

122.5 

75 

78.56 

37.41 

1.56 

2.0 

1.57 

10 

0.1 

193.75 

106.25 

43.9 

20.63 

4.41 

5.15 

4.65 

2 

1 

85 

40 

64.06 

26.45 

1.33 

1.51 

1.38 

2 

2 

78.75 

38.75 

65.47 

46.82 

1.20 

0.83 

1.05 

4 

1 

107.5 

72.5 

65.47 

42.09 

1.64 

1.72 

1.67 

4 

4 

92.5 

48.8 

64.53 

51.37 

1.43 

0.95 

1.22 

10 

1 

185 

78 

61.72 

31.8 

3.0 

2.45 

2.81 

10 

10 

122.5 

75 

78.56 

37.41 

1.56 

2.0 

1.57 

20 

1 

257 

112 

57.98 

29.93 

4.43 

3.74 

4.20 

20 

20 

150 

100 

55.25 

37.25 

2.71 

2.68 

2.70 
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Table  6:  Influence  of  interphase  region  on  composite  oxidation 


Di,  axial/ 

D  ox 

Di,trans/ 

D  ox 

taxial, 

laxial, 

ttrans,Hm 

ltrans, 

taxial/ttrans 

laxial/ltrans 

- 

- 

57.5 

23.75 

43.9 

20.6 

1.31 

1.15 

1 

l 

59.14 

22.18 

47.74 

23.87 

1.24 

0.93 

10 

l 

95.43 

32.26 

49.23 

22.38 

1.94 

1.44 
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Fig.  1.  Oxidation  measurement  procedures 


Fig.  2.  Photomicrograph  of  PMR-15  resin  showing  the  formation  of  oxidized  region 


center  (nonoxidized) 


9CUIUIICU 
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Fig.  3.  Evolution  of  oxidation  layer  and  transition  region  thickness  with  aging  time  in 

PMR-15  resin  at  288°C 


Fig.  4.  Oxidation  layer  thickness  for  ambient  and  elevated  pressure  aging  in  PMR-15 


Aging  time,  hrs 

Fig.  5.  Normalized  weight  loss  of  G30-500  fiber  and  PMR-15  resin 
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(b) 


Fig.  6.  (a)  Array  of  eight  stitched  photomicrographs  of  oxidized  specimen  cross  section, 


(b)  Enhanced  micrograph  of  the  specimen  cross  section 


Fig.  7.  Unidirectional  G30-500/PMR-15  composite  at  407,  1200  and  2092  hours  of 

oxidation  at  288°C 
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Fig.  8.  Comparison  of  oxidation  growth  in  axial  and  transverse  directions  of  composite 


Fig.  9.  Comparison  of  oxidation  growth  in  transverse  direction  of  composite  and  in  neat 

resin  PMR-15 
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Direction  of  oxidation  propagation 


oxidized/ 

degraded 

material 


active 

oxidation 

process 

zone 


unoxidized 

material 


Fig.  10.  Schematic  of  the  three-zones  in  thermo-oxidation. 


Fig.  11.  Prediction  of  oxidation  thickness  at  316°C  and  343°C  based  on  288°C  model 

parameters 
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Fig.  12.  Influence  of  sorption  on  the  oxidation  layer  thickness 


Aging  time,  hrs 


Fig.  13.  Prediction  of  oxidation  thickness  for  pressurized  aging  of  PMR-15 
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Tom  3D:  The  Thermo  Oxidation  Modeler  for  High  Temp.  PMfc  |fn~|^<] 


Fig.  14.  User  input  interface  for  TOM-3D 
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Surface  Exposed  to  02 


Surfaces  Exposed  to  02 


Fiber 

interphase 
domains  shown 
In  exploded  view 


(a)  (b) 

Fig.  15.  Representative  volume  elements  used  for  the  analyses 
(a)  Single  Fiber  RVE  (12.5  pm  x  12.5  pm  x  125  pm)  with  one  surface  exposed  to 
oxidation;  (b)  Ten  Fiber  RVE  (125  pm  x  12.5  pm  x  125  pm)  with  both  axial  and 
transverse  directions  exposed  to  oxidation. 
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(a)  (b)  (c) 

Fig.  16:  Effect  of  fiber  diffusivity  on  composite  oxidation.  Black  regions  are  completely 
oxidized  and  grey  regions  are  unoxidized.  The  transition  region  is  contoured  from  black 
to  grey  based  on  the  oxidation  state.  All  dimensions  shown  are  in  pm. 

(a)  Df  =  0.01  Dmox ,  (b)  Df  =  0.1  Dmox,  (c)  Df=  10  Dmox 
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10  hrs 


With  Interphase  (Df  =  0.1  Dmox;  Dijaxiai=  10  Dmox;  Dijtrans=  Draox)  :  taxia|/ttrans  =  1-935 


No  Interphase  (Df  =0.1  Dmox);  taxiai/ttrans  =  1-212 


Fig.  17.  Effect  of  a  highly  diffusive  fiber-matrix  interphase  region.  The  anisotropy  in 
the  oxidation  layer  can  be  clearly  seen. 
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